Shear Zone Refraction and Deflection in Layered Granular Materials 
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Refraction and deflection of shear zones in layered granular materials was studied experimentally 
and numerically. We show, that (i) according to a recent theoretical prediction [T. linger, Phys. 
Rev. Lett. 98, 018301 (2007)] shear zones refract in layered systems in analogy with light refraction, 
(ii) zone refraction obeys Snell's law known from geometric optics and (iii) under natural pressure 
conditions (i.e. in the presence of gravity) the zone can also be deflected by the interface so that 
the deformation of the high friction material is avoided. 



When granular materials deform under external stress 
the deformation is often localized into narrow regions. 
These shear zones [Il[ai3l[l[5l|6l[aia[9l[T0l[IIl[l2l 
[T3l HH |T5l [161 nil [IH] ^ct as internal slip surfaces be- 
tween solid-like blocks of the bulk. The formation of 
shear zones is a crucial deformation mechanism in fine 
powders, sand and soil (landslides). Geological faults are 
themselves large scale examples of shear zones. Here we 
study experimentally a recent theoretical prediction that 
shear zones behave in striking analogy with geometric 
optics [111 [13] , and compare our results with numerical 
simulations. We show that shear zones alter their orien- 
tations when crossing media boundaries similarly to light 
refraction but here the frictional properties of the mate- 
rials take the role of the optical refractive index. We find 
that the refraction phenomenon also exists in the pres- 
ence of gravity, i.e. under natural pressure conditions. In 
certain configurations we observe another effect, namely 
that shear zones can be deflected by the material inter- 
face. 

In the present experiments we use two materials with 
different frictional properties. One material is corun- 
dum which consists of angular grains therefore it has 
higher effective internal friction than the other mate- 
rial consisting of glass beads. To characterize the dif- 
ference in the internal frictions the angles of repose (Or) 
were determined by the method used in [19 and were 
Of'' = 21.9° for glass beads while l9^^^ = 33.2° for corun- 
dum which is somewhat higher than the typical value for 
sand (l9^^^ = 30.5°). The ratio tan l9^^7 tan l9^^^ = 1.63 
gives a reasonable contrast for the internal friction. The 
interface between the two materials (illustrated with dif- 
ferent colors in Figs.[TJa) and [ijb)) is a tilted plane. The 
inclination of this tilted plane with respect to horizontal 
defines the interface angle 0. Shearing was performed in 
a straight split bottom cell [H [6l [71 [10] with internal cross 
section of 4.2 cm x 4.5 cm. The shear cell included two 
60 cm long L shaped sliders (cell wall), one of which was 
slowly translated in the experiments according to the red 
(gray) arrows (with total displacement between 5 and 6 
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FIG. 1: (color online) Schematic illustration of (a) refraction 
and (b) deflection of a shear zone in layered granular systems 
according to our experimental setup, (c)-(d) Top view of the 
corresponding experiments after deformation. The measure- 
ment was done with colored and uncolored particles in order 
to visualize the displacement profile (illustrated by white ar- 
rows). (e)-(f) Shear strain at the top measured during the 
experiment. The dashed line shows the position of the split 
line at the bottom. 

cm), which creates the shear zone indicated by the red 
line on the sketch. The presence of the interface leads to 
refraction or deflection of the shear zone depending on 
the configuration. This is in contrast to the simple case 
when the cell is filled with one material and the shear 
zone is formed along a vertical plane [3l [6l [7| . 

The idea of refraction for shear zones was proposed 



earlier [12] and studied numerically in a special geometry 
with cylindrical walls and periodic boundary conditions. 
In that setup both ends of the shear zone were fixed and 
the system was sheared under constant pressure in zero 
gravity. These computer simulations showed that the 
zone changes its orientation at the interface and the ex- 
tent of refraction follows Snell's law of light refraction. 
With the present work we were able to reproduce the ef- 
fect of refraction experimentally, moreover we show that 
the phenomenon exists also under more realistic circum- 
stances. In our case the position of the shear zone is free 
at the top surface and gravity keeps the grains together 
leading to a natural pressure gradient. This geometry 
allowed us also to detect the effect of deflection where 
the interface modifies the position of the shear zone in 
such a way that the deformation of the high friction re- 
gion is completely avoided. Deflection was not observed 
in the original cylindrical cell [12 because the two ends 
of the zone were fixed on different sides of the interface. 
Therefore, the zone was forced to cross the interface and 
to enter the high friction material. 

The deformation in our case was measured as follows. 
For both materials we used two samples with different 
colors, therefore four regions appear at the top surface as 
it is seen in Figs.[TJc) and[TJd) . Two of these regions cor- 
respond to glass beads (grain size dp = 0.56 ±0.02 mm for 
the colored sample and dp = 0.48 ± 0.02 mm for the un- 
colored sample) while the other two correspond to corun- 
dum (grain size dp = 0.33 ±0.02 mm for the colored sam- 
ple and dp = 0.23 ± 0.02 mm for the uncolored sample). 
Using particles with different sizes allowed the separation 
of samples after each experiment, but have been shown 
not to affect the results [20j . This way not only the inter- 
face of the two materials is visualized (oriented vertically 
in Figs. [TJc) and [ijd)) but the displacement profile is 
also directly seen, as it is illustrated with white arrows 
in Figs. [TJc) and TJd). The deformation of the surface 
layer was recorded during translation by a video camera, 
movies recorded during the experiments can be seen at 
[20] . The displacement at the top surface as a function 
of the position y was determined by digital image analy- 
sis. The corresponding shear strains (see Figs. [TJe) and 
[TJf)) show that the shear zones are shifted away from the 
split at the base of the cell. They are not in the middle 
(at y = 0) where they would be in the absence of the 
interface i.e. for a homogeneous material. 

After each experiment the displacement profile in the 
bulk d{y^ z) was reconstructed by removing the top sur- 
face of the material carefully layer by layer [18 . We used 
a commercial vacuum cleaner with additional extension 
tubes to provide a slow flow rate. The procedure takes 
several hours for one experiment. The displacement pro- 
files are presented in Figs. |2|a) and |3|a) for refraction 
and deflection, respectively together with sample images 
taken during the excavation process. The gradient of 
the displacement provides the local shear strain 7 inside 
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FIG. 2: (color online) (a) The reconstructed displacement 
profiles in the bulk for refraction for = 74° taken by care- 
fully removing the material layer by layer (sample images 
taken during excavation are shown on the left hand side). 
Color (grayscale) corresponds to the deformation d{y,z). 
Movies demonstrating the excavation process can be seen at 
[20 . (b) The corresponding shear strain distribution in the 
y — z plane. The colors show the dimensionless shear strain 
7 as defined in the text. The tilted blue (dark gray) line 
indicates the position of the interface. 



shear strain 7 as the function of position (Figs. [2|b) and 
|3jb)) the structure of the shear zone becomes visible. In 
Fig. |2|b) the zone, starting from the bottom, takes a 
short path towards the interface and by reaching the low 
friction region it changes direction abruptly and heads 
straight to the top. In Fig. |3|b) where the zone starts 
in the low friction material it deflects to avoid the high 
friction region even if it takes a much longer path. 

The behavior can be well understood qualitatively 
based on a simple variational model of shear zones 
[T2I [16]. This model regards the shear zone as an in- 
finitely thin slip surface and states that the shape of the 
slip surface is chosen according to minimum energy dissi- 
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FIG. 3: (color online) (a) The reconstructed displacement 
profiles in the bulk for deflection for = 59° taken by care- 
fully removing the material layer by layer (sample images 
taken during excavation are shown). Movies demonstrating 
the excavation process can be seen at [20]. (b) The distribu- 
tion of the shear strain 7 in the y — z plane. The tilted blue 
(dark gray) line indicates the position of the interface. 



pation. For our straight geometry this is reduced to the 
condition that the minimization of the following integral 
provides the shape of the zone: 



jj.p dl 



(1) 



where the local effective friction coefficient ja times the 
local pressure p is integrated along the path that is taken 
by the shear zone in the ?/-z-plane section. Thus not 
the length of the path but the length weighted by jap 
is minimized by the shear zone. This is quite similar to 
Fermat's principle of optics where the length weighted by 
the refraction index is minimized by light beams. This 
analogy leads to the idea that similar refraction effects 
were expected in these two distant fields of physics. 
We performed experiments for numerous orientations 
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FIG. 4: The ratio sin a/ sin/3 {a and p being the angles of 
incidence) for zone refraction as a function of the interface an- 
gle (j) with respect to horizontal obtained from (a) experiments 
and (b) numerical simulations. For the case of simulations the 
open data points were taken in the absence of gravity (i.e. ho- 
mogeneous pressure) . The horizontal dashed lines correspond 
to the ratio of the effective frictions /icor/Mgia = 1-63 deter- 
mined from the angles of repose for the two materials used. 
The inset of (b) illustrates the angles a, f3 and 0. 



of the interface to test whether zone refraction obeys 
Snell's law known from optics, which reads sin a/ sin/3 = 
Mcor/ Mgia in the present case. Here a and /3 are the angles 
of incidence (see inset of Fig. ^h) for illustration) and 
the relative index of refraction is replaced by the ratio 
of the effective frictions /icor//^gia [12] • Fig. [4]^ a) shows 
sin a/ sin/3 as a function of the interface angle (j) with 
respect to horizontal and the value of /icor/^gia = 1-63 
(dashed line) that was determined experimentally based 
on the two angles of repose [19 . As it is seen the data 
are independent of (j) (within errors) for a wide range of (j) 
and scatter around the expected value 1.63 in agreement 
with Snell's law. 

In order to better compare the predictions of the above 
concept with our experimental data we performed com- 
puter simulations based on the fluctuating narrow band 
model [17 . This method follows the concept of minimum 
dissipation but also takes the fact into account that the 
optimal slip surface has to be found in a random medium. 
Technically, in the simulation we search for an optimal 
path in a two dimensional underlying lattice which rep- 
resents the cross section of the shear cell. Statistical 
fluctuations are introduced by random strength of the 
lattice-bonds. Depending on the random realization of 
the medium different slip surfaces are obtained. Taking 
ensemble average over random realizations one arrives at 
a continuous displacement profile instead of a thin slip 
surface. To incorporate the ratio of the effective fric- 
tions /icor/ Mgia we use the same ratio for the average bond 
strength between the low and the high friction regions. 
There is one fit parameter in the simulation, the extent 
of fluctuations, which is chosen so that the width of the 
shear zone at the free surface approximately matches the 
experimental value. More details on this method can be 
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found in [TT. 

Numerical simulations (Fig.|4jb)) confirm with higher 
accuracy that sin a/ sin/3 matches the expected value of 
/icor/^gia = 1.63. Simulations also revealed, that this 
is true for both the absence and the presence of gravity 
where the pressure was taken to be constant and pro- 
portional to the depth, respectively |3j. In gravity the 
shear zone is slightly curved which is a direct consequence 
of the pressure gradient and makes the measurement of 
refraction more difficult. Under homogeneous pressure 
where the curvature is absent, the accuracy of data is 
much higher (see open circles in Fig. [4jb)). 

For a more detailed quantitative comparison two sim- 
ulations are shown in Fig. [5] which correspond to the ex- 
perimental configurations presented in Figs. |2] and [3j The 
shear strain distributions obtained by the simulation are 




FIG. 5: (color online) Distributions of the shear strain 7 ob- 
tained in numerical simulations for (a) refraction and (b) de- 
flection based on the fluctuating narrow band model. The 
dashed lines show the center of the experimentally achieved 
shear zones, for comparison. 

in very good agreement with the experiments. The shape 
of the shear zones and other details are nicely reproduced, 
e.g. the largest values of the shear strain. In the case 
of refraction the largest strain is obtained at the bottom 
and right above the interface near to the refraction point. 



while for defiection the high shear strain is located at the 
bottom and at the top of the material. This latter effect 
leads to a narrower shear zone at the top in our second 
configuration (defiection) than in the first one (refrac- 
tion), which is also nicely detected in the experiments 
- compare the width of the zone in Figs. [TJe) andjljf). 
Comparing the center positions of the shear zones at the 
top gives also nice agreement. We obtain y = 0.85 cm 
in the simulation and 0.81 cm experimentally in case of 
refraction, while for defiection the positions are —1.11 cm 
and —1.15 cm in the simulation and experiment, respec- 
tively. 

In this work we studied shear zones in granular sys- 
tems where two different materials are layered on top of 
each other. On one hand we provided an experimental 
evidence of the phenomenon of zone refraction, on the 
other hand we found a new effect where the shear zone 
is defiected by the interface. For the case of zone refrac- 
tion a close analogy can be drawn with light refraction 
in optics. Our experimental and numerical data show, 
that the well known Snell's law is valid for this granu- 
lar system. The presented experiments and simulations 
clearly demonstrate that both effects - refraction and de- 
fiection - alter significantly the position of the shear zone, 
thus they may play important role in various industrial 
and civil engineering applications. Moreover, we expect 
that our findings may be relevant in geological processes 
where the materials involved (soil, rock) have naturally 
layered structures. 

Further experiments in layered materials are on the 
way to explore the bulk deformation with an indepen- 
dent method (MRI imaging). We recently learned about 
a parallel experimental study on shear zone refraction 
obtained in a vertical cylindrical cell [21]. 
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